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Singular and Half-Quantum Vortices and Associated Majorana Particles in 
Superfluid ^He-A between Parallel Plates 

Takuto Kawakami*, Yasumasa TSUTSUMI, and Kazushige Machida 

Department of physics, Okayama University, Okayama 700-8530 

(Received December 16, 2009) 

Motivated by a recent experiment on superfluid "^He-A confined in narrow parallel plates 
using a rotating cryostat, we explore possible vortices stable under magnetic field applied to 
arbitrary angle relative to the plates in order to seek vortices which can accommodate the 
Majorana zero mode in the core. After proving that the singular vortex with the unit winding 
number provides the Majorana mode in the spinful situation, we establish the phase diagram 
in the plane; the rotation frequency Q vs system size R by finding possible order parameter 
textures within the Ginzburg-Landau framework. We also analyze the stability for a single and 
a pair of half-quantum vortices, which possesses the Majorana mode in its core. It is concluded 
from the above mention that the Majorana zero mode can be found in the present on- going 
experimental setting at ISSP, Univ. Tokyo. 

KEYWORDS: ^He, A-phase, chiral p-waye pairing, superfluid, parallel plate geometry, vortex, texture 



1. Introduction 

There has been much attention on topological orders 
and the associated Fermionic quasi-particles with low 
energies. This is particularly true when the quasi- 
particles are Majorana^^ character, namely its creation 
operator 7''' is equal to its annihilation operator 7; 7''' = 
7. This unusual and intriguing Fermion, quite differ- 
ent from the usual Dirac particle, is thought to be use- 
ful for fault-tolerant quantum computations because it 
obeys non-Abelian statistics^^ and its existence is pro- 
tected topologically to avoid decoherence. These situ- 
ations are ideal for quantum computation.^^ The can- 
didate systems, which support the Majorana particle, 
are quite rare; chiral spinless Px ± ipy superconductors, 
p-wave Feshbach resonanced superfluid, and the frac- 
tional quantum Hall state with the 5/2 filling. The former 
superconductors have not been identified yet in nature. 
It has often been argued that Sr2Ru04 may be a prime 
candidate, ^^^^^ but strong doubt has been cast on this 
possibility of Sr2Ru04 of its triplet pairing. ^^^^^ Note 
that the first discovered triplet superconductor UPta is 
an /-wave pairing, not chiral p-wave.^^^^^ It is proven 
theoretically that the half-quantum vortex (HQV) in the 
chiral superconductors with the p-wave pairing possesses 
the Majorana particle with zero energy localized at the 
vortex core.^^'^"*^^ 

Superfluid ^He-A phase is characterized by a chiral 
p-wave pairing. There is no doubt on this identifica- 
tion.^^' In fact, Volovik and Mineev^^^ are the first 
to point out the possibility to the realization of HQVs 
in 1976. Since then, there have been several general ar- 
guments on the stability of a HQV in connection with 
^He-A phase. ^^^^^ However, there are no serious calcu- 
lations which consider it in a realistic situation in su- 
perfluid ^He-A phase on how to stabilize it and on what 
boundary conditions are needed for it. 

Recently, Yamashita, et al?"^^ have performed an ex- 
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periment intended to observe HQVs in superfluid ^He-A 
in parallel plate geometry. The superfluid is confined in 
a cylindrical region with the radius = 1.5 mm and the 
height 12.5 jiia sandwiched by parallel plates. A mag- 
netic field H = 26.7 mT {\\z) is applied perpendicular to 
the parallel plates under pressure P=3.05 MPa. Since the 
gap 12.5 jam between plates is narrow compared to the 
dipole coherence length ~ 10/im, the /-vector, which 
signifies the direction of orbital angular momentum of 
Cooper pairs, is always perpendicular to the plates. Also 
the d-vector is confined within the plane by an applied 
field H \\ z because the dipole magnetic field Hd ~ 2.0 
22,23) tends to align the (i- vector perpen- 

dicular to the field direction. They investigate to seek 
out various parameter spaces, such as temperature T, or 
the rotation speed 1] up to 11 = 6.28 rad/sec by using 
the rotating cryostat in ISSP, Univ. Tokyo, capable for 
the maximum rotation speed ~12 rad/sec, but there is 
no evidence for the existence of the HQV.^^^ 

The aims in this paper are to investigate the possible 
vortex structures which can accommodate the Majorana 
particle in the core under the above realistic experimental 
situations at ISSP for superfluid ^He-A phase confined 
in the parallel plates. The candidate vortex structures in 
this situation are either the singular vortex (SV) with 
odd winging number or HQVs mentioned above. Thus 
after examining the sufficient condition for the Majorana 
zero energy particle to exist in the SV, we determine 
the phase diagram in the system size with the radius 
R and the external rotation frequency Q under applied 
fields with arbitrary angle relative to the plates. Note 
that the external field is necessary for performing NMR 
detection. 

The arrangement of the paper is as follows: In §2 we 
examine the possible order parameters and its spatial 
structures, or textures. By utilizing general symmetry 
properties of the superfluid ^He-A phase, we demonstrate 
that the SV with odd winding number even in the spinful 
situation can accommodate the Majorana zero mode in 
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Table I. The combination of phase windings of each OP compo- 
nent for the A-phase texture (AT), the singular vortex (SV) and 
the half-quantum vortex (HQV). 



the core. In §3 we investigate the stable vortices and tex- 
tures within the Ginzburg-Landau (GL) formahsm when 
the magnetic field is applied to an arbitrary angle rela- 
tive to the plates, in order to find the stable SV region 
for various system parameters, external rotation the 
system size R and the field orientation Oh- In §4 the sta- 
bility problem of the single HQV and a pair of HQVs 
are analyzed within the same framework. We devote to 
summary and conclusions in the final section. 
A part of the §4 is published in ref. 30. 

2. General Considerations 

2. 1 Order parameters and textures 

Generally, in p-wave superfiuids, the order parameter 
(OP) is described by 

A||(r,p) A||(r,p) 
The matrix element is described by 



A(r,p) 



(1) 



.{r)p,. (2) 

where a is spin index | or | , j9± = ^(pa^ ± ipy)/^/2^ 
and p is the unit vector in the momentum space. In the 
parallel plate geometry, the momentum pz component is 
suppressed. Here each component Acrcr',±(^) {o- =t, i) 
can have its own phase winding whose winding num- 
ber is denoted by Wcra' ,±' There are three possible tex- 
tures at rest and under rotation as shown in Table I: 
(K;crcr',+7 '^crcr',-) = (O^^); A-phasc texture (AT), (1,3); 
SV, and (^||,+ , ^||,+, _, _) (0,1,2,3); HQV. 
These textures and vortices are allowed in axisymmetric 
situation (see §3.2). We investigate the detailed configu- 
ration of these textures, the relative energetics, and the 
associated Majorana quasi-particle in the last two tex- 
tures. 

2.2 Majorana hound state -symmetry considerations - 

In this subsection, we examine the existence of the 
Majorana particle in the SV and the HQV in a most 
general situation. By assuming that the configurations of 
the OP are the bulk A-phase, the SV texture is described 
as 

^aa'ir.p) = A^^/,+(r)j9+ 

= d^^/(r)A+(r)exp(i0)p+, (3) 

where (j) is the azimuthal angle of the center-of-mass co- 
ordinate r. 



/-vector 




^-vector \^ ^ \^ 



Fig. 1. (Color online) Schematic diagram of the parallel plate sys- 
tem. The external magnetic field H tilts by Oh from the z-axis. 
Od is the polar angle of the vector. The external field tends to 
align the ci-vector perpendicular to H and the dipole interaction 
tends to aligned the vector parallel to the direction z. Then 
< < 7r/2 - Oh. 



and vector d = {dx^dy^dz) is real. Under this assump- 
tion, the direction of vector is determined by the dipole 
interaction and the interaction with the external field. 
The polar and azimuthal angle of direction of d-vector 
are determined as 



I (bulk) 



tan 



sin 26 H 



{Hd/HY-cos2eH 



0. 



(4) 
(5) 



The angle Od and Oh are defined as shown in Fig. 1. 
The angle (j)^ is the azimuthal angle of the direction of 

vector. In other words, when the direction of the spin 
quantization axis is in the angle of {Q^^^^^^ +7r/2, ^^^^^^^), 
du (^) = throughout the system. 

When A||(ri,r2) = A||(r,p) = 0, p-wave mean field 
Hamiltonian is written as 



H = y dridr2^\rx) 



/Cu 



*(r2), 



(6) 



where 



-A*,(ri,r2) 



2m3 



Aaa(^l,^2) 

-i7o"*(n,r2) 



^(^l) - Ma 



-lii [xi- yi^— 

oyi oxi 



S{ri - r2), 



where ma, and /icr are the mass of the ^He atom, the 
single particle potential, and the chemical potential of 
the particle whose spin a. We carry out the Bogoliubov 
transformation to quasi particles whose creation and an- 
nihilation operators 
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j dr [vyi{r)il;^{r) + Vy2{r)il^i{r) 

+Uyi {r)^lj\ (r ) + (r')'0 j (r) . 
Obtained Bogoliubov-de Gennes (BdG) equation is 



(7) 



(8) 



%(ri,r2) 



2*1. (^i) 





El 



(9) 



wiiere = (73 This (73 is the Pauli matrix. We can 
reduce the BdG equation (9) to four eigenvalue equa- 
tions, 

Uyi{r2) 
Vvi{ri) 



j c^r2/C||(ri,r2) 
j dr2JCii{ri,r2) 



(10) 



(11) 



These eigenvalue equations are solved numerically.^^^ If 
Acrcr(^,p) has odd phase winding, we can obtain the 
zero energy state {E^ = 0) and {uj^^i{r),Vjy^i{r)) = 
('?;* •(r),7i*^^(r)) so that r]j,a = That is, the Bogoli- 
ubov quasi-particles exhibit the Majorana character. 

In the SV texture all components of the OP have the 
odd winding number. Thus the SV has the Majorana 
quasi-particle with the external field toward any direc- 
tion. In the HQV texture, either A||(r,p) or A||(r,p) 
has the odd phase winding. The component that has no 
phase windings does not have the lower excitation. That 
is, the HQV has the lower excitation that is induced by 
the odd phase winding and these excitations contain the 
Majorana quasi-particle. Consequently, both the SV and 
the HQV have the Majorana quasi-particle at these vor- 
tex cores. 

However, in the above analysis, we do not consider the 
influence of the vortex core and the edge of system. In 
these situation, the configuration of the OP is different 
from the bulk A-phase due to variance of the OP, mix- 
ture of minor components, and so on. Furthermore, the 
Majorana quasi-particle exists at the vortex core and the 
edge of system. Then we have to clarify OP textures by 
more realistic and serious calculation using GL theory. 



3. Ginzburg-Landau Functional and Phase Dia- 
gram 

3. 1 Ginzburg-Landau functional 

The GL free-energy functional invariant under gauge 
transformation, spin, and orbital space rotations is well 
established 22,23,27,32 37) given by a standard form 

f total ~ fgrad ~l~ fbulk ~l~ fdipole ~l~ f field') (1^) 

fsrad = K [{d:A;j){diA^j) + {d:A;^){djA^i) 

+ {dtA*MdjA^,)] , 



fbulk 



'^P^^lii^ti^l-^j^i^j + (^"^^lii^Zj^lJij^i^i 



fdipole — 9d{^ijiiji^vu ~\~ ^ij,iy^iyiJ. q ^/uz^^Atz^) 5 



f field — QmHfiA^j^HiyAi^i^ 

where = x,y,z and di = \/i — {2ims/h){ft x r)i 
{ft II z) with 7713 being the mass of the ^He atom. The 
coefficients of quadratic terms 

aaa' = c^o(l - T/Tc + sgn(a + a')^T/T^), (13) 

where = A^(0)/3. The Zeeman effect splits the tran- 
sition temperature Tc| of the jH) pair and of the 
lii) pair.^s) The split width AT = (Tct - Tc|)/2 in eq. 
(13). However, we consider the region that AT <C Tc — T 
so-called weak field region until §4. The coefficient of 
fourth order terms are obtained in ref. 33 appro- 

priate for the experiment at P = 3.05 MPa. We use 
the gradient coupling constant with weak coupling limit 
K = 7C{3)N{0){hvFf /240{7rkBTc). The g/^^ and the 
gm are the coupling constant of the dipole interaction 
and the interaction with external field respectively. 



Qd 



Mo 
40 



jhN{0) In 



1.1339 X 0.45Ti. 



(14) 



(15) 



48[{l^F-)7rkBT,]' 

The transition temperature under no magnetic field Tc, 
the density of states A^(0), Fermi velocity vp^ the perme- 
ability of vacuum /io, the gyromagnetic ratio 7, Landau 
parameter Tq , and Fermi temperature Tp are given by 
experiments.^^' In order to carry out realistic calcu- 
lations taking account of the configuration of the vortex 
core and the edge of the system, our formulation includes 
strong coupling corrections in the bulk terms where the 
A-phase is stable over the B-phase. Thus the gradient 
terms are within the weak coupling limit {psp/ps = 1) 
because there is no established method to properly take 
into account the Fermi liquid correction for the OP form 
generalized beyond the A-phase (see for standard method 
ref. 23). We minimize the GL functional eq. (12) where 
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the external field is applied to arbitrary direction relative 
to the plates and absolute value. First, we consider two 
limits where the direction of the external magnetic field 
is parallel and perpendicular to the plates. 

3.2 Parallel field H\\ 

In the situation where the external field applies par- 
allel to the plates, both the dipole interaction and the 
external field suppresses the components of the OP 
A||(r) and A||(r). In the bulk A-phase system, the two 
components A^i^±{r) survive. Since the phase winding 
of the OP components increases the kinetic energy of 
Cooper pairs, it is energetically favorable at rest that 
large component of the OP has small winding num- 
ber. Furthermore, in the cylindrical symmetric system, 
w^l^j^ = 1^11- —2 due to orbital coupling effect. Under 
low rotations the stable texture is the AT whose phase 
winding {w^i^j^.,w^i^-) = (0,2) and the dominant com- 
ponent of the OP is A^i^-^{r). 

Under the anti-clock wise rotation the component that 
has the positive phase winding gains the rotating energy 
through a term — i^iiilA^iip in the GL functional eq. 
(12). As the rotating speed Q is higher, the free energy of 
the SV whose dominant component of the OP has wind- 
ing number w = 1 becomes lower than that of the AT. 
According to the above analysis, the candidates of stable 
SVs are (^u,+.^n,-) = (-1.1) and (1,3). The (-1,1) 
texture has smaller winding number than the (1,3) tex- 
ture, from the kinetic energy view point the (—1,1) tex- 
ture seems to be more stable than the(l,3). However, 
under rotation the minor component A||_(r) of (1,3) 
texture gains more rotating energy than the minor com- 
ponent A||+(r) of (—1,1) texture. Under rotation the 
stable texture is determined by competition of those two 
factors. The results of our calculation show that when 
the rotating speed is large enough to stabilize the SV 
not AT, (1, 3) texture is always more stable than (—1, 1) 
texture. 

In Fig. 2, the cross section of the amplitude of the 
OP component of the AT and the SV when the direc- 
tion of the external field is parallel to the plates and 
the quantization axis of spin is chosen perpendicular to 
the plates. Choosing the quantization axis parallel to the 
plates same as the direction of the field, we can obtain 
the SV texture whose OP component A||(r) = ex- 
actly. The SV texture has Majorana quasi-particle state 
at the vortex core. 

The phase diagram of textures in the space consist- 
ing of the rotating speed Q and system size R is shown 
in Fig. 3 in a fixed temperature T/Tc = 0.95. We calcu- 
late in various temperatures and absolute values of exter- 
nal field, showing that the temperature and the absolute 
value of external field change Qc only slightly, so that Qc 
is determined by almost only R. We obtain the critical 
rotating speed Qc quantitatively as a function of system 
size R smaller than 100 jam. Therefore, we can extrap- 
olate Qc = 0.06 rad/sec in R = 1.5 mm, that is the 
system size of the sample using the experiment in ISSP. 
This rotating speed can be well controlled by the present 
experimental technique. 



0.5 

I 



AT 



SV 



f WU+ = 1 



' w^i- = 3 . 



r[|im] 



5 



r [|am] 



Fig. 2. (Color online) Order parameter amplitude nor- 
malized by TvkBTc. Left (right) figure shows the cross section of 
the OP components along the radial direction r of the AT (SV) 
texture forR = 5 /im, T/Tc = 0.95, H = 10 mT, Oh = 7t/2. We 
choose the spin quantization axis perpendicular to the plates; 
the z direction. The spin component Aero- vanishes. 




10 100 1000 

System size R [|im] 

Fig. 3. System size R dependence of the critical angular velocity 
Qc from the AT to the SV T/Tc = 0.95, H \\ plates {Oh = 7t/2). 
An extrapolated value of flc = 0.06 rad/sec is found at R = 1.5 
mm. The inset shows the energy differences Af between the two 
textures as a function of for R = 100 fim. 



3. 3 Perpendicular field H± 

In this subsection we consider the situation where the 
external field apply parallel to the plates. When the ex- 
ternal field larger than the dipole field H^, the (i- vectors 
lie parallel to the plates. The result of numerical calcu- 
lation shows that A||(r,p) = exactly throughout the 
system as shown in Fig. 4. These SV textures have the 
Majorana quasi-particle state at the vortex core. When 
the external field is smaller than the dipole field Hd, in 
the bulk A-phase the (i- vector aligns perpendicular to the 
plates. The OP texture is same as in the case that the 
external field is parallel to the plates shown in Fig. 2. 
Choosing the quantization axis parallel to the plates, we 
can obtain the OP component A^i{r) = exactly. How- 
ever, when vector || H taking account of Zeeman effect 
in BdG equation (9), the symmetry of the equation ar- 
rows the solution {uj,^i{r),Vj,^i{r)) = (^^,2(^)5 ^^,2(^)) • 
Thus the core bound state of the SV cannot satisfy 
Majorana condition {uj,^i{r),Vj,^i{r)) = {vl^i{r),ul^^{r)) . 
Therefore only in the case where the external field is 
larger than the dipole field Hd^ the stabilized SV has the 
Majorana quasi-particle. 

The phase diagram of textures are same as that in 
§3.2 since in the result of our calculation, the external 
magnetic field change critical rotating speed Vtc of the 
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Fig. 4. (Color online) Order parameter amplitude \A^^ 
malized by TrksTc forR = 5 /im, T/Tc = 0.95, H = 
6h = 0. Left (right) figure shows the cross section of the OP 
components along the radial direction r of the AT (SV) texture. 
We choose the spin quantization axis perpendicular to the plates; 
z direction. The spin component A||(r,p) = exactly. 
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Fig. 5. (Color online) Order parameter amplitude |A^cr',±| nor- 
malized by T^kBTc for = 5 /xm, T/Tc = 0.95, H = 2 mT, 
6h = 7r/18. Left (right) figure shows the cross section of the OP 
components along the radial direction r of the AT (SV) texture. 
We choose the spin quantization axis perpendicular to the plates; 
z direction. All spin components A^^f(r,p) are non- vanishing. 
When we choose the spin quantization axis to be the direction 
e = ef""^^^ + 7r/2 and (/> = 0, < Aaa,±. 



texture transition from the AT to the SV only slightly as 
shown in §3.4. 

In the case that H ± plates where the external field 
is larger than the dipole field Hd^ the vector is per- 
pendicular to the /-vector. Then the dipole energy of the 
bulk system can be neglected. Thus this case is favorable 
for realization of the HQV texture. We discuss the HQV 
texture in §4 in detail. 

3.4 Arbitrary oriented field 

In the bulk A-phase system, the direction of the d- 
vector is determined by the competition between the 
dipole interaction favoring that the (i- vector is parallel to 
/-vector and the external field favoring that the vector 
is perpendicular to the H. The direction of the vector 
is given by eqs. (4) and (5). When the quantization axis 
is perpendicular to the plates, all components Acrcr'(^,p) 
are non- vanishing as shown in Fig. 5. 

We consider the energetics of these textures by calcu- 
lating the stable textures under various external mag- 
netic fields; = 1, 2, 3, and 10 mT, for various polar 
angles 6h = nn/lS n = 1,2,---. The result of the calcu- 
lation in the system R = 5 /im which is shown in Fig. 6. 
Under these external fields, Qc coincides with each other 
for the precision in second order after the decimal point. 
We see that the external magnetic field scarcely changes 
the phase diagram shown in Fig. 3. Their reasons are 
as follows. The critical rotation Qc is determined by the 
competition of the kinetic energy and the energy gain 



Fig. 6. (Color online) The critical angular velocity flc depen- 
dence of Oh from the AT to the SV for T/Tc = 0.95, R = 5 
/im, if = 1, 2, 3, and 10 mT. Under the various direction and 
absolute value of the external fields, flc coincides with each other 
for the precision in second order after the decimal point. 



from angular momentum, namely, the gradient term. On 
the other hands the infiuence of the interaction with the 
external field is no more than 10~^ order of the gradient 
energy. 

We discuss the possibility that the Majorana quasi- 
particle exists in these textures. As shown in §2.2, the 
SV has the Majorana quasi-particle when their OP com- 
ponent A||(r,p) = for appropriate quantization axis 
throughout the system. Thus we turn the quantization 
axis to be their polar and azimuthal angle (0,0) ^ 
^Q^huik) ^ (j)^^^^^^) in order to change the expression 
to that suppressing the component of the OP A||(r,p). 
Then we obtain changed OP textures as shown in Fig. 7. 
In these textures, the components of the OP A^i^±{r) 7^ 
0, but finite amplitude with the order 10~^ against the 
dominant component Acrcr±(^)- These amplitudes vary 
spatially, and their phase cannot be specified by the ac- 
curate A-phase OP. In addition, when the vector is not 
perpendicular to the magnetic field we should note 
the infiuence of Zeeman effect for the Majorana condi- 
tion. Therefore, we cannot conclude whether or not those 
SVs have the Majorana quasi-particle. In order to clar- 
ify these points, we have to solve spinful BdG equation 
include off-diagonal terms of eq. (9). These remain as a 
future problem. 

We notice these induced components. In Fig. 8, we 
show the maximum value of amplitude of induced com- 
ponent A^i^j^ir) as a function of 9h in various exter- 
nal field values. We show that when the direction of the 
external field are perpendicular to the plates {Oh = 0) 
and parallel to the plates {6h = 7r/2), the component 
A||^+(r) = exactly. Moreover, when the external field 
comparable to the dipole field Hd 2 mT), the induced 
component is enhanced maximally under such a magni- 
tude of the field. 

The system size dependence of them is shown in Fig. 
9. When the system size becomes large, the amplitude 
of the component A^i^-^{r) tend to be large. Therefore, 
even if the large sample such as that used in experiment 
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Fig. 7. (Color online) Order parameter amplitude (a) |Acrcr,+ |, (b) |Acrcr,-|, (c) and (d) normalized by tt/c^Tc for 

R = 5 /im, T/Tc = 0.95, H = 2 mT, Oh = 7r/18. We choose the spin quantization axis = 0^^^^^^ -\- cj) = 0. The spin components 
A||(r,p) / but have the amplitude with the order 10~^ of the amplitude of Aaa(r,p) and modulate spatially. 




9^[rad] 




System size R [|im] 



Fig. 8. (Color online) The maximum value of as a func- 

tion of Oh for H = 2, 3, and 27 mT, R = 5 /im, T/Tc = 
0.95. When Oh = or 7r/2, the amplitude of the component 
A||(r,p) = exactly When the external field is comparable 
to dipole field (^ 2) mT, the amplitude of these components is 
enhanced. 



Fig. 9. The maximum value of as a function of R for 

H = 1,2, and 3 mT, Oh = 7t/6 and tt/S rad, /im, T/Tc = 0.95. 
When the system size becomes larger, the amplitude of these 
components becomes larger. In the system using experiment, 
R = 1.5 mm, are non-vanishing. 



R = 1.5 mm, these components are finite. 

We suggest ttiat ttie reason wiiy ttie induced com- 
ponent is non-zero is due to the existence of induced 
Acra' component. Taking account of the A^a',- compo- 
nent, the direction of vector determined by minimizing 
f dipole + f field has the Component parallel to the spin 
quantization axis except for the case that H \\ plates 
and H ± plates. In addition, these components varies 
spatially according to spatial modulation of A^^i _ com- 
ponent. 

4. Half-Quantum Vortex 

J^.l Single half-quantum vortex 

When the direction of the external field is perpen- 
dicular to the plates, it is a favorable situation for the 
HQV to realize. Then in this section, we consider pos- 
sible stable HQV texture and their energetics. First, we 
consider the case that there is a single HQV in the sys- 
tem. The external magnetic field splits Tc of the |tT) 
and III) pairs by the Zeeman effect. The splitting is 
defined as AT = (T^t -Tc|)/2. When AT > (AT < 0), 
the amplitude of the OP A^^ (A||) becomes larger than 
A|| (A||). The splitting AT is generally much smaller 



than the amplitude of the OP in the so-called weak field 
region :^ 10 mT. In this region, the texture of the 
HQV are shown in Fig. 10. The rotation speed Q^ci (S^c2) 
at which the free energy of the HQV (SV) and the AT 
(HQV) are intersect. In Fig. 11, we show that the HQV 
never becomes the absolute stable texture. The reasons 
are the strong coupling effect in the bulk fourth terms of 
GL functional eq. (12). In the weak coupling limit, under 
some critical rotation, the AT, the SV, and the HQV de- 
generate each other. The strong coupling effect stabilize 
the A phase, then the structure of the central region of 
the AT are more favorable than the others. Therefore, 
the critical rotation l^ci are larger than ^c2- 

When the temperature becomes larger, the amplitude 
of the OP becomes smaller. Then the split AT cannot 
be neglected. When the split AT > 0, the |tT) pair be- 
comes dominant. Thus, the HQV texture and the en- 
ergy difference A/ in Fig. 11 become close to the AT. 
In this process, l^ci become smaller than ^c2- Further- 
more, when the temperature becomes larger, the stable 
region of the HQV becomes larger as shown in Fig. 12. 
Consequently, the phase diagram is shown in Fig. 13. 
We can extrapolate that the stability region of the HQV 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 7 



0.5 


AT 




sv 




|Tt> 




|U> 




^w^^^ = 


\ 




^ mi+ = 1 


\ 




wtt- = 2 













r [|im] 


5 


r [|im] 


5 



□ An+ 



Fig. 10. (Color online) Order parameter amplitude of (a) A^^^± 
and (b) A^^^± components normalized by tt/c^Tc for = 5 /xm, 
T/Tc = 0.95, if = 10 mT, Oh = 0. We choose the spin quantiza- 
tion axis perpendicular to the plates; z direction. At the vortex 
core, component is only non- vanishing, so that the HQV 

is A 1 -core vortex. 
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Fig. 11. (Color online) Free-energy comparison for the AT, the 
HQV, the SV, the two HQVs system (2HQV), and the two SVs 
system (2SV) as a function of Q for R = 10 /im and T/Tc = 0.97 
at the weak field region. The Af is the relative free energy to 
the AT. In the weak field region, Qc2 < ^ci- 
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Fig. 12. (Color online) Stability region of the HQV in fl versus 
T/Tc (R = 10 /im and AT/Tc = 0.05). AiT denotes the Ai 
phase texture where only |tT} pair exist. The stability region of 
the HQV becomes wide at the vicinity of the transition temper- 
ature from the A-phase to the Ai -phase. 



is 0.05 < Q < 0.06 rad/sec in the sample whose size 
R = 1.5 mm. We also notice that by changing the radius 
R of the system one can control the width of the stability 
region. For example, in R = 100 /im, the stability region 
of the HQV is 7 < 1] < 8 rad/sec. 

4.2 A pair of half- quantum vortices 

We examine the case that there are more than two 
half-quantum vortices in the system. There are two pos- 
sibilities to enter the two HQVs in the system. Namely, 
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Fig. 13. Stabihty region of HQV sandwiched between Qci and 
nc2 as a function of R for T/Tc = 0.97 and AT/Tc = 0.05 
(H ~ 100 mT). The critical rotation flci = 0.05 rad/sec is the 
extrapolated value for R = 1.5 mm. Inset shows the free-energy 
comparison for R = 20 /im, displaying the successive transitions 
from the AT to the HQV at Qd and from the HQV to the SV 
at Qc2- 



the one case is that there are the HQVs whose wind- 
ing number , = (0,1,0,3) and 
(1,0,3,0) in the system, the other case is that there are 
two (0,l,2,3)-HQVs. We calculate these two cases and 
conclude that the former is not solution of our calcula- 
tion. The center of phase winding of the component of 
the OP does not coincide with the center of the system, 
which is not the energetically advantageous form for the 
angular momentum. Furthermore, there is no repulsion 
between HQVs because we assume that the gradient en- 
ergy are the weak coupling limit form. 

We investigate the case that there are two (0,1,2,3)- 
HQVs in the system. In the Fig. 14 we show the texture 
of the two HQVs system (2HQV). In this texture, the 
component that has the phase winding are the same spin 
state. Thus the distance between them is determined by 
the competition of the kinetic energy loss and the energy 
gain from the angular momentum. However, we conclude 
these are just metastable texture from our calculation. 

First of all, we consider the case that AT < where 
the III) pair is dominant. In this case, the texture of 
the 2HQV is close to two SVs system (2SV). From the 
analogy with the single HQV case, it can be guessed that 
critical rotation ^(at^2HQV) where the free energy of the 
2HQV cross to that of the AT is larger than the critical 
rotation ^(2HQV^2SV) where the free energy of the 2SV 
crosses to that of the AT. Namely, the 2HQV cannot be 
stable in this course. 

On the other hands, we consider the case that AT > 
where the |||) pair is dominant. In this case the texture 
of the 2HQV are energetically close to the AT. From 
the analogy with the single HQV case, it is possible that 



^(AT^2HQV) 



< n 



(2HQV^2SV) 



. However, in these regions. 



the SV is more stable than the AT. From our calculation, 
we show that the critical rotation speed ^(hqv^2HQV) 
where the free energy of the 2HQV crosses to that of 
the HQV cannot smaller than the critical rotation where 
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Fig. 14. (Color online) Order parameter amplitude (a) A^^^± and 
(b) components normalized by tt/c^Tc for the met ast able 

2HQV system, R = 5 fim, T/Tc = 0.95, Oh = 0. We choose the 
spin quantization axis perpendicular to the plates; z direction. 
There are two vortices whose winding number 1 in the \ space. 
The sense of the phase winding and its number are shown. The 
low figures are cross sections of each component across the vortex 
cores. 
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either exactly parallel or perpendicular to the plates, the 
singular vortex state becomes stable above the critical 
rotation speed Vtc which is evaluated. Thus the Majo- 
rana quasi-particles can be observed in those situation. 
The on-going experiment using the rotating cryostat at 
ISSP, Univ. Tokyo is particularly suited for this task be- 
cause Vtc — 0.06 rad/sec is within well controlled rotation 
speed. Moreover since singular vortex carries the Majo- 
rana quasi-particle localized at the core, we can design 
the braiding experiment by manipulating two or more 
Majorana quasi-particles. 

Our calculations on the relative stability of the relevant 
textures are based on the standard Ginzburg-Landau 
functional which is firmly established through the cross 
checking between experiments and theories over several 
decades. Thus we believe that our results should be reli- 
able not only qualitatively, but also quantitatively. 

There are a few question remained to be clarified; (1) 
When the field is applied to exactly neither 7^ 0, or 
7^ 7r/2, there appears the unwanted OP component 
whose magnitude is an order of 10~^ compared with 
the main component. At present we do not know whether 
or not this component is really harmful for the existence 
of the Majorana quasi-particle for the singular vortex. 
(2) We do not know how to take into account the so- 
called strong coupling effect in the gradient term in eq. 
(12) for general order parameter state. This remains a 
future problem. 

In conclusion, independent of the above two reserva- 
tions (1) and (2) we emphasize the fact that the parallel 
plate geometry whose gap is an order of 10 jiia for super- 
fiuid ^He-A phase is one of the best systems for detecting 
the Majorana quasi-particle. 
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Fig. 15. Critical angular velocity ^(hqv^2HQV) from the HQV 
to the 2HQV for = 5 ^uv and AT/Tc=0.05. Even if the tem- 
perature becomes higher, ^(hqv^2HQV) does not become lower 
than r2c2- 



the free energy of the SV crosses to that of the HQV 
as shown in Fig. 15. Namely the 2HQV cannot also be 
stable in this course. 

5. Summary an conclusions 

We have studied the relative stability of the three kinds 
of the order parameter textures; the A-phase texture, the 
singular vortex, and the half-quantum vortex for super- 
fiuid ^He-A phase confined by narrow parallel plates with 
the radius R. After showing in §2 that the spinful singu- 
lar vortex can accommodate the Majorana quasi-particle 
localized at the core, we examine the condition that the 
Majorana quasi-particle exists, namely, the stability of 
the singular vortex and also half-quantum vortex which 
is known to accommodate the Majorana quasi-particle. 

We have found that under the external field applied to 
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